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Abstract. We study M-tensors and various properties of M-tensors are given. Spe- 
cially, we show that the smallest real eigenvalue of M-tensor is positive corresponding to a 
nonnegative eigenvector. We propose an algorithm to find the smallest positive eigenvalue 
and then apply the property to study the positive definiteness of a multivariate form. 
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1. Introduction 

M-matrices [H [121 EE] are known to have many applications in modeling dynamic sys- 
tems in economics, ecology, and engineering [23j. Various properties of M-matrices were 
used in establishing stability results for dynamic systems in general [23l[2l]. Especially, it 
follows from the famous Perron- Frobenius theorem for nonnegative matrices that for any 
M-matrix A, all real eigenvalues of A are positive fi2[ 126] . The theory of M-matrices 
has many applications in many fields such as computational mathematics, mathematical 
physics, mathematical economics, wireless communications, etc. [H |26]. Since avoidance 
conditions were linked to a stabihty type of result via Lyapunov-type functions, the theory 
of M-matrices was also used to certify avoidance conditions |25]. On the other hand, testing 
positive definiteness of a multivariate form is an important problem in the stability study 
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of nonlinear autonomous systems via Lyapunovs direct method in automatic control [T6] . 
Researchers in automatic control studied the conditions of such positive definiteness inten- 
sively [H 121 El [H]. For n > 3 and m > 4, this problem is a hard problem in mathematics. 
There are only a few methods to solve the problem [21 |3l [16]. In practice, when n > 3 and 
m > 4, these methods are computationally expensive. Note that an eigenvalue method was 
proposed in |16] to solve the problem. In [13], a method for computing the largest eigenvalue 
of an irreducible nonnegative tensor was applied to test the positive definiteness of a class of 
of multivariate forms. Motivated by these applications and methods, we study higher-order 
M-tensors in this paper. We show that all real eigenvalues of an M-tensor are positive. 
Hence, testing the positive definiteness of an even-order multivariate form is equivalent to 
testing a tensor is an M-tensor. 

A tensor can be regarded as a higher-order generalization of a matrix, which takes the 
form 

^=iAi-im), Ai^...i^ e R, l<ii,...,im<n. 

Such a multi-array A is said to be an m-order n- dimensional square real tensor with n™' 
entries Ai^...i^. In this regard, a vector is a first-order tensor and a matrix is a second- 
order tensor. Tensors of order more than two are called higher-order tensors. M-tensor we 
defined in this paper is ultimately related to the nonnegative tensor. It is a higher order 
generalization of the so-called M-matrix. 

Nonnegative tensors, arising from multilinear pagerank [13], spectral hypergraph theory 
[1], and higher-order Markov chains [15], etc., form a singularly important class of tensors 
and have attracted more and more attention since they share some intrinsic properties 
with those of the nonnegative matrices. One of those properties is the Perron- Frobenius 
theorem on eigenvalues. In [5], Chang, Pearson, and Zhang generalized the Perron- Frobenius 
theorem from nonnegative matrices to irreducible nonnegative tensors. Later, Yang and 
Yang [21] generalized the weak Perron-Frobenius theorem to general nonnegative tensors, 
and proved that the spectral radius is an eigenvalue of a nonnegative tensor. Numerical 
methods for finding the spectral radius of nonnegative tensors are subsequently proposed. 
Ng, Qi, and Zhou [15j provided an iterative method to find the largest eigenvalue of an 
irreducible nonnegative tensor. The Ng-Qi-Zhou method is efficient but it is not always 
convergent for irreducible nonnegative tensors. Zhang and Qi [29] established an explicit 
linear convergence rate of the Ng-Qi-Zhou method for essentially positive tensors [19]. Liu, 
Zhou and Ibrahim [T3] modified the Ng-Qi-Zhou method such that the modified algorithm is 
always convergent for finding the largest eigenvalue of an irreducible nonnegative tensor. The 
linear convergence rate of the algorithm was established in [30] for weakly positive tensors. 
Chang, Pearson and Zhang [^ introduced primitive tensors and established convergence 
of the Ng-Qi-Zhou method for primitive tensors. Friedland, Gaubert and Han [9] pointed 
out that the Perron-Frobenius theorem for nonnegative tensors has a very close link with 
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the Perron-Frobenius theorem for homogeneous monotone maps, initiated by Nussbaum 
[T7t [18] and further studied by Gaubert and Gunawardena [10] . Friendland, Gaubert and 
Han [9] gave a weaker definition for irreducible nonnegative tensors, and estabhshed the 
Perron-Frobenius theorem in this context. 

By using the eigenvalue theory of nonnegative tensors, we give some properties of M- 
tensors. We prove that the smallest eigenvalue of an M-tensor is positive with a nonnegative 
eigenvector. Let ^ be a tensor with nonpositive off-diagonal entries. Two necessary and 
sufficient conditions for A as an M-tensor are given. Specially, we prove that A is an M- 
tensor if and only if its all real eigenvalues are positive. We also give a sufficient condition 
for a tensor to be an M-tensor, which is easily verified. Finally, we propose an algorithm 
for computing the smallest real eigenvalue. The proposed algorithm is always convergent 
for any M-tensor. Furthermore, we link M-tensors with a class of multivariate forms and 
then apply the proposed method to study the positive definiteness of a multivariate form. 
It should be pointed out that the class of multivariate forms studied in [H] is a special case 
of our model. We do not need the assumption that the diagonal entries are positive. 

This paper is organized as follows. In Section 2, we will recall some preliminary results. 
We will introduce M-tensors and characterize some basic properties of M-tensors in Section 
3. In Section 4, we will propose an iterative algorithm for finding the smallest real eigenvalue 
of an M-tensor, and some numerical results are reported. In Section 5, we will present an 
algorithm to for testing positive definiteness of a class of multivariate forms. We conclude 
the paper with some remarks in Section 6. 

2. Preliminaries 

We start this section with some fundamental notions and properties on tensors. An m- 
order n-dimensional tensor A is called nonnegative (or, respectively, positive) if Ai-- i„i > 
(or, respectively, Ai^...i^ > 0). The tensor A is called symmetric if its entries Ai-^...i^ are 
invariant under any permutation of their indices {ii - ■ ■ im} [20] • The m-order t?,- dimensional 
unit tensor, denoted by X, is the tensor whose entries are with Sij^...i^ = 1 if and only if 

ii = ■ ■ ■ = im and otherwise zero. A tensor A is called reducible, if there exists a nonempty 
proper index subset Jc {1,2,. such that 

= 0, Vzi G /, Vz2, . . . , V ^ 
Otherwise, we say A is irreducible. 

Analogous with that of matrices, the theory of eigenvalues and eigenvectors is one of the 
fundamental and essential components in tensor analysis. Wide range of practical applica- 
tions can be found in [131 1211 122]. Compared with that of matrices, eigenvalue problems for 
higher-order tensors are nonlinear due to their multilinear structure. 
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Definition 2.1 Let C be the complex field. A pair {X,x) G Cx (C"\{0}) is called an 
eigenvalue- eigenvector pair of A, if they satisfy: 

Ax""-^ = Xx^"'-^\ (1) 

where n- dimensional column vectors Ax"^^^ and x'™"^' are defined as 

•^•^ •= I ^ ] ■^ii2---im-^i2 ' ' ' -^im I and X^ ' .= )l<j<n' 

\i2,...,jm = l / l<j<„ 

respectively. 

This definition was introduced in |6l [131 ISO]- We call p(^) the spectral radius of tensor A if 

p{A) = max{|A| : A is an eigenvalue of A}, 

where |A| denotes the modulus of A. An immediate consequence on the spectral radius 
follows directly from Corollary 3 in [20] . 

Lemma 2.1 Let A be an m-order n-dimensional tensor. Suppose that B = a{A + bl), 
where a and b are two real numbers. Then ^ is an eigenvalue of B if and only if = a{\ + b) 
and A is an eigenvalue of A. In this case, they have the same eigenvectors. 

Clearly, from this lemma, we have p{B) < \a\ {p{A) + \b\). 

Nice properties such as the Perron- Frobenius theorem for eigenvalues of nonnegative 
square tensors have been established in [S]. 

Theorem 2.1 If A is an irreducible nonnegative tensor of order m and dimension n, then 
there exist Aq > and Xq > 0, Xq E such that 

Moreover, if A is an eigenvalue with a nonnegative eigenvector, then A = Aq. If X is an 
eigenvalue of A, then |A| < Aq. 

Clearly, by this theorem, Aq is the largest eigenvalue of A. 

Yang and Yang [2H] asserted that the spectral radius of a nonnegative tensor is an 
eigenvalue. In the following we state the results of [28t Theorem 2.3 and Lemma 5.8] for 
reference. 
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Theorem 2.2 Assume that A is a nonnegative tensor of order m dimension n, then p{A) is 
an eigenvalue of A with a nonzero nonnegative eigenvector. Moreover, for any x > 0, x G i?" 
we have 

{Ax'"~^)i ^ ... {Ax"'-^)i 



— < P(^) < max 



mm , ^ f., ^ ^ ..^^v^ 1 



The foUowing continuity of the spectral radius was given in the proof of [281 Theorem 2.3]. 

Lemma 2.2 Let A be a nonnegative tensor of order m and dimension n, and e > be a 
sufficiently small number. If Ae = A + S where S denotes the tensor with every entry being 
e, then 

\imp{Ae) = p{A). 

£->-0 

Now we state an iterative algorithm for calculating the largest eigenvalue of a nonnegative 
tensor A, which is proposed in [TH [15] based on Theorems 12.11 and 12.21 

Algorithm 1 

Step 0. Choose a;(^) > 0, x^^'> G R'' and p>0. LetB = A + pi, and set k:=l. 
Step 1. Compute 



, m—1 



Ak = mm^(fe)^p— — — 



Afc = max (k) 



x) 



Step 2. If \k = Xk> 'l^hen let X = \k and stop. Otherwise, compute 



X 



{k- 



replace k by k + 1 and go to Step 1. 
Step 3. Output X — p, the largest eigenvalue of A. 



It is proved in [H] that Algorithm [T] is always convergent for irreducible nonnegative 
tensors; see the following theorem. 



5 



Theorem 2.3 Suppose A be an irreducible nonnegative tensor. Let B = A + pi, where 
p > 0. Assume that A is the largest eigenvalue of B. Then, AlgorithmU\ produces a value of 
X in a finite number of steps, or generates two sequences {A^} and {Xk} which converge to 
A. Furthermore, X — p is the largest eigenvalue of A. 

Note that for any nonnegative tensor .4, ^ + £^ is an irreducible nonnegative tensor. 
Therefore, if we set i3 = ^ + pX + S in Algorithm [H then by Lemma [2.21 we can prove that 
the modified algorithm is also convergent in the similar way. 

3. M-tensors 

We now extend the notion of M-matrices to higher-order tensors and introduce the 
definition of an M-tensor. 

Definition 3.1 Let A be an m-order and n-dimensional tensor. A is called an M-tensor 
if there exist a nonnegative tensor B and a real number c > p{B), where p{B) is the spectral 
radius of B, such that 

A = cl-B. 

Clearly, when m = 2, if ^ is an M-tensor then A is an M- matrix. 

Note that the off-diagonal entries of an M-tensor are nonpositive. Denote Z the set of 
m-order and n-dimensional real tensors whose off-diagonal entries are nonpositive. 

Note that eigenvalues defined in Definition 12. II exist for an M-tensor [20j. We will show 
that every eigenvalue of an M-tensor has a positive real part, and hence all of its real 
eigenvalues are positive. 

Theorem 3.1 Let A be an M-tensor and denote by t{A) the minimal value of the real 
part of all eigenvalues of A. Then t{A) > is an eigenvalue of A with a nonnegative 
eigenvector. Moreover, there exist a nonnegative tensor B and a real number c > p{B) such 
that t{A) = c — p{B). If A is irreducible, then t{A) is the unique eigenvalue with a positive 
eigenvector. 

Proof. Since A is an M-tensor, by Definition 13.11 there exist a nonnegative tensor B and 
a real number c > p{B) such that 

A = cl-B. 

Let A be an eigenvalue of A and ReA be the real part of A. By Lemma 12. 1^ c — A is an 
eigenvalue of B. Since p{B) is the spectral radius of B, 

p{B) > |c - A| > c - ReA > p{B) - ReA, 
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which imphes that Re A > 0, and hence 

p{B) > max {c - ReA} = c - min {ReA} = c - t(A). (2) 

On the other hand, Theorem 12.21 shows that p{B) is an eigenvalue of B. By Lemma 12. 1^ 
c — p{B) is a real eigenvalue of A, and hence 

c-p{B)>r{A), 

which, together with ([2]), implies 

T{A) = c-p{B). 

Therefore, t{A) is an eigenvalue of A. Since p{B) has a nonnegative eigenvector x*, x* is 
also an eigenvector of t{A). Moreover, if A is irreducible, then B is also irreducible. Hence 
p{B) is the unique eigenvalue of B with a positive eigenvector. Thus, we complete the proof. ■ 



We state a few of the alternative necessary and sufficient conditions of a tensor in Z to 
be an M-tensor. 

Theorem 3.2 A tensor in Z is an M-tensor if and only if any of its eigenvalues has a 
positive real part. 

Proof. Let ^ G Z, and suppose that every eigenvalue of A has a positive real part. Let 
a = max {Ai i}. Then B = aX — ^ is nonnegative. By Lemma 12.11 and Theorem 12.21 

l<j<n 

a — p{B) is a real eigenvalue of A. It follows that a — p{B) > 0. Thus ^ = al — i3 is an 
M-tensor. 

Now, let ^ = cX — i3 be an M-tensor. Then, B is nonnegative and c > p{B). Let A be an 
eigenvalue of A. Let ReA be the real part of A. Suppose ReA < 0. Since A is an eigenvalue 
of A, there exists a nonzero vector a; G C" such that 

which yields 

Bx""-' = (c- A)xl'"-^l 

It follows from Definition 1.1 that c — A is an eigenvalue of B. Since c > p{B) and ReA < 0, 
we have 

|c - A| > c - ReA > c> p{B), 
which contradicts the maximality of p{B). We complete the proof. ■ 



Another equivalent condition for an M-tensor is presented as follows, which is more easily 
certified than the one in Theorem 13.21 
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Theorem 3.3 A tensor A E 7^ is an M-tensor if and only if all of its real eigenvalues are 
positive. 

Proof. The necessity of the condition follows directly from Theorem 13.21 Suppose now that 
all real eigenvalues of ^ G Z are positive. Let B = aX — A where a = max {Ai i}. Then B 

l<i<n 

is nonnegative and hence it follows from Theorem 12.21 that p{B) is an eigenvalue of B. By 
Lemma \2A\ a — p{B) is a real eigenvalue of A, it must be positive. That is, a > p{B), and 
therefore A is an M-tensor. ■ 



Clearly, we easily obtain the following results from the above theorems. 
Corollary 3.1 If A E 7^ is an M-tensor, then max{y4j...j} > 0. 

l<i<n 

Proof. Let a = maxjAj j} and B = aX — A. Then i5 > and hence p{B) is an eigenvalue 

l<i<n 

of B. So, a — p{B) is a real eigenvalue of A. Since A is an M-tensor, by Theorem 13.31 
a — p{B) > 0, which implies a > p{B) > 0. ■ 

Note that. If the tensor A E 7 then A can be deposed into the form A = XX — B with 
A > max{y4j.. .(} and B>Q. Thus, we have the following sufficient and necessary condition 

l<j<n 

for a tensor being an M-tensor, which is a generalization of Theorem 7 in |27]. 

Corollary 3.2 Let A E 7. Decompose the tensor A into the form A = XX — B, where 
A > max {A. ..?.}• Then A is an M-tensor if and only if X > p{B). 

l<i<n 

Proof. Clearly, B > 0. By Theorem 12. 2[ p{B is an eigenvalue of B. By Lemma |2TT| X — p{B) 
is a real eigenvalue of A. If A be an M-tensor then we have, by Theorem 13. 3[ A — p{B) > 0, 
i.e., A > p{B). If A > p{B) then, by Definition 13. A is an M-tensor. ■ 

Note that this result provides us an easy method for determining whether a tensor A 
is an M-tensor. We only need to compute the spectral radius of the tensor XX — A, where 
A > max{Ai i}. 

l<i<n 

Finally, we give a sufficient condition for a tensor to be an M-tensor, which is easily 
verified. First, we state a definition which is a generalization from matrices to tensors 

[HIES]. 

Definition 3.2 Let A be an m-order and n-dimensional tensor. A is diagonally dominant 
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n 

^ ] |^M2...im I — I^M---i|) i = l,2,...,n, (3) 

where the symbol 5^2. ..i™ = defined as the entry of the unit tensor I. A is strictly 
diagonally dominant if the strict inequality holds in ^ for all i. A is irreducibly diagonally 
dominant if A is irreducible, diagonally dominant, and the strict inequality in ^ holds for 
at least one i. 



Theorem 3.4 // a tensor A E 7^ is strictly or irreducibly diagonally dominant with all 
nonnegative diagonal entries, then A is an M-tensor. 

Proof. Let A be an eigenvalue of A with a nonzero eigenvector x. Denote 

\Xi\ = max . 

Then 

n 

i2v)*m = l 

which imphes that 



n 



^^2 vi^m — ^]^ii2 . ..im — ^ 



Hence, the diagonal dominance of A implies that 

|ReA - A, J <\X- A, J < \Ai,J. (5) 
Since > for j = 1,2, ... ,n, ([5]) yields 

ReA - A,„i > -A,„i. (6) 

Suppose that A is strictly diagonally dominant. Then the strict inequality holds in ([3]) 
for all j, so the strict inequality holds in (Q. This yields ReA > 0, i.e., any eigenvalues of A 
has a positive real part. By Theorem 13. 2 [ ^ G Z is an M-tensor. 

Suppose now that A is irreducibly diagonally dominant. Define 
J = {I : \xi\ = max \xi\, \xi\ > |xj|for some i}. 

l<i<n 

If J = 0, then (jl]) and the diagonal dominance of A imply that for i = 1, 2, . . . , n, 

n 

I \ _ /I .. .\ < \^ M.. . \ < \ -I 

i2,.--,*m = l;(5ii,...i™ =0 
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Let 

n 

\^kk...k\ > ^ ] I^A:i2...iml 

*2v>«m = l;'5fci2.-.im=0 

for some k. We have 

|ReA — Ak...k\ < |A — Ak...k\ < \A.k...k\ = A.k...k 

which imphes that Re A > 0. 

If J 7^ 0, then the irreducibihty of A imphes that there exist / G J and i2, . . ■ ,im ^ J 
such thatv4/i2...j^ 7^ 0. Hence (jl]) yields 

" I I I I ^ 

\^-Aii,j\< \ l^ii2...j™lTn rn~ ^ I^«i2...iml < 

i2,...,«m = l;<5li2...im=0 «2 , ■ ■ ■ ,*m = 1 ;0;i2 . . .im =0 

which imphes that ReA > 0. By Theorem 13.21 ^ G Z is an M-tensor. ■ 



4. An Algorithm 

Theorem 13.31 shows that a tensor ^ G Z is an M-tensor if and only if the smallest real 
eigenvalue of A is positive. In this section we propose an algorithm to determine whether 
or not a tensor with nonpositive off-diagonal entries is an M-tensor. 

Lemma 4.1 Let A be an m-order and n- dimensional tensor. Define 

Lj\^= mm{Aii,„i-Ci}, [/^ = max {Ai...* + C'i}, (7) 

l<i<n l<i<n 

where 

n 

Ci = ^ ^ 1^222. ..im 1 1 ^ 1, 2, . . . , 71. 

i2vi*m = l,5ii2---»m=0 

Then and Ujs^ are the lower and upper hounds of real eigenvalues of A, respectively. 
Proof. Let A be a real eigenvalue of A with an eigenvector a; 7^ 0. That is, 

n 

Aii2...im^i2' ' ' ^im = ^' i=l,2,...,n. (8) 

«2,---i«m = l 

Let [xfcl = maxi<i<„ \xi\. Then ([8]) implies that 

n I ' I I ■ I 

|A — v4fcfc...A:| < |^fc»2..iml'r~V • • • I '"i ^ Ck, 
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which yields Akk...k — Ck < X < Akk...k + C*fc- This shows < A < U^. m 

For a tensor ^ G Z, we define a tensor C as 

C = Uj(I - A, (9) 

where Uj^ is defined in ([7]). Clearly, the tensor C is a nonnegative tensor. By Lemma WA\ and 
Theorem 12.11 Ujs, — p{C) is the smallest eigenvalue of A. By Theorem 13.31 and Definition 13. 1[ 
if t/4 — p{C) is positive then A is an M-tensor. Based on this observation, in the following, 
we establish an algorithm for computing the smallest eigenvalue of a tensor A E 'Z. If the 
smallest eigenvalue is positive then A is an M-tensor. Otherwise, A is not an M-tensor. 

Based on the above discussion, we propose the following algorithm to determine whether 
or not a tensor A with nonpositive off-diagonal entries is an M-tensor. 

Algorithm 2 

Step 0. Compute the upper bound of real eigenvalues, U_a by the formula (0j and let C = 
UaI - A. 

Step 1. By using AlgorithmUl compute X, the largest eigenvalue of C. 

Step 2. Let p = f/4 — A. If fi > then A is an M-tensor. Otherwise, A is not an M-tensor. 

In order to show the viability of Algorithm [21 we used Matlab Version 7.7(R2008b) to 
test it on some tensors with nonpositive off-diagonal entries which are randomly generated 
by the following procedure. 

Procedure 1. 

1. Give {m,n, Ad), where n and m are the dimension and the order of the randomly 
generated tensor, respectively, and > 0. 

2. Randomly generate an m-order n-dimensional tensor V such that all elements of V are 

in the interval (0,1). 

3. Let A = where A-i = + A-i, i = 1,2, otherwise, Ai-i™ = 

Our numerical results are reported in Table[T]In this table, n and m specify the dimension 
and the order of the randomly generated tensor, respectively. A^ is a parameter in Procedure 



11 



1. Given {m,n, Ad), we generate 100 tensors and determine whether or not they are M- 
tensors by Algorithm [2l In the Yes column we show the number of tensors which are 
M-tensors. In the No column, we give the number of tensors which are not M-tensors. 
CPU(s) denotes the average computer time in seconds used for Algorithm El The results 
reported in Table [T] show that Algorithm [2] performs well for these tensors. 

5. An Application: The positive definiteness of a multivariate form 

In this section we apply the proposed algorithm for testing the positive definiteness of a 
class of multivariate forms. 

An mth degree homogeneous polynomial form of n variables f{x), where x G -R", can be 
denoted as 

n 

/(x) .= ^ ^ Ai-^^i^ i^Xij^Xi^-'-Xi^. (10) 

ii,i2,...,im=l 

When m is even, f{x) is called positive definite if 

f{x) > 0, Vx G R", x^O. (11) 

Testing positive definiteness of a multivariate form is an important problem in the stability 
study of nonhnear autonomous systems [21 |3l [16]. It is proved in [22] that f{x) is positive 
definite if and only if its corresponding tensor A = (Aj^ j^) is symmetric and all of real 
eigenvalues are positive. That is, testing positive definiteness of a class of even-order mul- 
tivariate forms is equivalent to determine whether or not the even-order symmetric tensor 
is an M-tensor. It should be pointed out that the positive definiteness of a special class 
of multivariate forms was studied in [14]. It can be regarded as our special case. Some 
numerical results can be referred [H]. We omit them here. 

6. Conclusions 

We have defined and studied M-tensors. An M-tensor is the generation of an M-matrix. 
Many important characterizations of M-matrices has been extended to M-tensors. We have 
proposed two sufficient and necessary conditions, a necessary condition, and a sufficient 
condition of an M-tensor. We have shown that the smallest eigenvalue of an M-tensor is 
positive, and then proposed an algorithm to determine whether or not a tensor with nonpos- 
itive off-diagonal entries is an M-tensor. Finally, we link an M-tensor with a multivariate 
form and apply the algorithm to judge the positive definiteness of the multivariate form. 
Numerical results are reported. 
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m n All 


Yes No CPU(s) 


3 10 5 
3 10 10 
3 10 100 
3 10 1000 


100 0.0025 
100 0.0030 
100 0.0028 
100 0.0028 


3 20 5 
3 20 10 
3 20 100 
3 20 1000 


100 0.0094 
100 0.0102 
100 0.0088 
100 0.0084 


3 30 5 
3 30 10 
3 30 100 
3 30 1000 


100 0.0181 
100 0.0187 
100 0.0186 
100 0.0195 


3 40 5 
3 40 10 
3 40 100 
3 40 1000 


100 0.0352 
100 0.0358 
100 0.0362 
100 0.0355 


3 50 5 
3 50 10 
3 50 100 
3 50 1000 


100 0.0619 
100 0.0619 
100 0.0598 
100 0.0692 


4 10 5 
4 10 10 
4 10 100 
4 10 1000 


100 0.0592 
100 0.0603 
100 0.0611 
100 0.0620 


4 20 5 
4 20 10 
4 20 100 
4 20 1000 


100 0.2945 
100 0.3097 
100 0.3187 
100 0.3116 


4 30 5 
4 30 10 
4 30 100 
4 30 1000 


100 1.3233 
100 1.3125 
100 1.3170 
100 1.3453 


4 40 5 
4 40 10 
4 40 100 
4 40 1000 


100 6.5375 
100 6.5358 
100 6.4925 
100 6.5520 


4 50 5 
4 50 10 
4 50 100 
4 50 1000 


100 15.2086 
100 15.1844 
100 15.2102 
100 15.2039 



Table 1: Output of Algorithm [2] for randomly generated tensors. 
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There are some questions are still in study. For example, whether the condition "there 
exists X e R'^,x > such that Ax^'^ > 0" is a necessary and sufficient condition for a 
tensor ^ e Z to be an M-tensor? 
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